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Abstract. We report on 49 fast-mode forward shocks propagating inside coronal mass 
ejections (CMEs) as measured by Wind and ACE at 1 AU from 1997 to 2006. Compared 
to typical CME-driven shocks, these shocks propagate in different upstream conditions, 
where the median upstream Alfven speed is 85 kms“^, the proton /3 = 0.08 and the 
magnetic field strength is 8 nT. These shocks are fast with a median speed of 590 kms“^ 
but weak with a median Alfvenic Mach number of 1.9. They typically compress the mag¬ 
netic field and density by a factor of 2-3. The most extreme upstream conditions found 
were a fast magnetosonic speed of 230 kms“^, a plasma {3 of 0.02, upstream solar wind 
speed of 740 kms“^ and density of 0.5 cm“^. Nineteen of these complex events were as¬ 
sociated with an intense geomagnetic storm (peak Dst under —100 nT) within 12 hours 
of the shock detection at Wind, and fifteen were associated with a drop of the storm¬ 
time Dst index of more than 50 nT between 3 and 9 hours after shock detection. We 
also compare them to a sample of 45 shocks propagating in more typical upstream con¬ 
ditions. We show the average property of these shocks through a superposed epoch anal¬ 
ysis, and we present some analytical considerations regarding the compression ratios of 
shocks in low (3 regimes. As most of these shocks are measured in the back half of a CME, 
we conclude that about half the shocks may not remain fast-mode shocks as they prop¬ 
agate through an entire CME due to the large upstream and magnetosonic speeds. 


1. INTRODUCTION 

East magnetosonic forward shocks are measured at Earth 
preceding fast coronal mass ejections (CMEs) and, more 
rarely, propagating into the slower stream of corotating in¬ 
teraction regions (CIRs). Many shocks are also observed 
without clear drivers, which probably reflects the fact that 
the angular extent of shocks is larger than that of the pu¬ 
tative, associated CME. Hereafter, we use CME to refer to 
coronal mass ejections, whether when they are observed re¬ 
motely, in the corona and in the heliosphere, or when they 
are measured in situ. During solar cycle 23, 250-300 shocks 
were measured by ACE and Wind, just upstream of Earth 
[e.g. Wang et al., 2010; Janvier et al, 2014]. Interplan¬ 
etary shocks are followed by a sudden storm commence¬ 
ment (SSC) or a sudden impulse (SI), where the magnetic 
field strength of low-latitude ground magnetometers suffers 
a sudden increase, typically associated with a compression 
of the magnetosphere by the solar wind dynamic pressure, 
and an increase in the Dst index. The sheath region behind 
a shock may be geo-effective; in fact, about 15% of intense 
geomagnetic storms (as defined by Dstpeak < —100 nT) in 
solar cycle 23 were caused by CME sheaths [e.g., see Zhang 
et al, 2007, and references therein]. Typical slow solar wind 
conditions at 1 AU are Alfvenic speed of ~ 50 kms“^, den¬ 
sity of 5 cm“® and solar wind speed of 400 kms“^ [see, 
for example Schwenn, 2006; Ebert et al, 2009]. Most inter¬ 
planetary shocks at 1 AU are relatively weak with a Mach 
number in the range of 2-4 and compression ratio of the 
order of 2 [e.g., see Berdichevsky et al, 2000; Oh et al., 
2007] for CME-driven shocks. In a survey of 105 quasi¬ 
perpendicular shocks, Richardson and Cane [2010b] (there¬ 
after, RC2010) found that about 30% of such shocks were 
propagating through or closely following a CME, represent¬ 
ing a fraction 2-3 times more important than what can be 
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expected from the proportion of CMEs in the solar wind 
[which is about 10-15%, as found in Richardson et ah, 2000; 
Zhao et ah, 2009]. 

The fact that shocks can propagate through CMEs has 
been known for a long time, with detections reported in 
Ivanov [1982] and Burlaga et al. [1987] for complex distur¬ 
bances measured in August 1972 and April 1979, respec¬ 
tively. Wang et al. [2003c] described in detail two occur¬ 
rences during solar cycle 23 associated with intense geomag¬ 
netic storms. The authors found evidence that the com¬ 
pressed magnetic held of the CME contributed to the en¬ 
hanced geo-effectiveness [Wang et al, 2003b]. Simulations 
of shocks propagating inside a flux rope or magnetic cloud 
have been performed by Vandas et al. [1997] and Xiong et al. 
[2006], among others. Lugaz et al. [2005] performed a numer¬ 
ical simulation of the interaction of two CMEs using three- 
dimensional magneto-hydrodynamic modeling, and studied 
in detail the propagation of the shock driven by the second 
CME inside the first magnetic ejecta, including the variation 
of the upstream conditions, the shock speed and compres¬ 
sion ratio (their Figure 2). They found that, in the main 
part of the cloud, where P ~ 0.1, the Alfven speed was 
~ 250 kms“^ and the sound speed was ~ 70 kms~^. Con¬ 
sequently, the shock compression ratio was found to decrease 
from about 3 before the CME and to reach 1.5 ±0.2 inside 
the CME. Based on the low compression ratio and the fact 
that magnetic ejecta are commonly regions of low density, it 
is unlikely that shocks inside CMEs may be imaged remotely 
by wide-angle white-light imagers such as the His onboard 
STEREO. In addition, remote-sensing does not give enough 
information about the upstream and downstream conditions 
to determine the shock properties and characteristics. A re¬ 
cent example of a shock propagating inside a CME as mea¬ 
sured in situ was reported for the current solar cycle by Liu 
et al. [2014]. 

Here, we expand on the work from RC2010 by study¬ 
ing in greater detail the occurences of shocks propagating 
through a CME during solar cycle 23 (1997-2008). We do 
so by combining the list of CMEs measured by Wind of 
Richardson and Cane [2010a] with the list of shocks mea¬ 
sured by Wind as identified and analyzed by J. Kasper and 
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Figure 1. Wind observations of the shock inside a CME 
on 6 August 1998. The slow CME starts around 13:00UT 
on 08/05 and the shock passes at 07:16UT on 08/06 
(31.26 h). The panels show from top to bottom, the pro¬ 
ton density (in red the alpha-to-proton number density 
ratio), the proton velocity and temperature (in red, the 
expected proton temperature), the total magnetic held 
and its components in GSM coordinates, the magnetic 
held clock angle, the Alfven speed, the Kan-Lee coupling 
function, the proton /3 and the Sym-H index. The CME 
boundaries as dehned by Richardson and Cane [2010a] 
are shown with black lines, the shock with a red dashed 
line. 

M. Stevens at the Center for Astrophysics (CfA) and the 
list of shocks measured by ACE as identihed at the Univer¬ 
sity of New Hampshire (UNH) by researchers and graduate 
students there [see, for example Paulson et aL, 2012]. Our 
goal in this endeavor is threefold: (1) to hnally dispose of the 
misconception that shocks do not propagate through CMEs, 
(2) to study the additional geo-effectiveness brought about 
by the shock’s compression of the CMEs, and (3) to study 
shocks in a regime which is unique to CMEs. To build upon 
this last point, we note that studying shocks inside CMEs 
give us a unique opportunity to have direct measurements 
of fast shocks in a magnetically dominated plasma (where 
the proton /3 is typically less than 0.1) and where the Alfven 
speed may reach values of several hundreds of kms“^. In 
section 2, we present various examples of shocks propagating 
through CMEs. In section 3, we present an overview of the 
properties of the shocks, the upstream conditions they prop¬ 
agate into, as well as the geomagnetic effects attributed to 
the compression of the CME. We also present a superposed 
epoch analysis of the 49 shocks. We discuss analytical con¬ 
siderations regarding the compression ratio of MHD shocks 


in low-/? plasma in section 4, and discuss our results and 
conclude in section 5. 

2. Examples of Shocks Inside CMEs 

2.1. Typical Case: August 6, 1998 shock 

We first present a typical case, chosen so as to be as close 
as possible to the median from our sample. Eigure 1 shows 
from top to bottom, the proton density, speed and tem¬ 
perature, the magnetic field strength, its components and 
clock angle (6), the Alfven speed, the Kan - Lee coupling 
function [Kan and Lee, 1979], the proton /?, and the Sym-H 
index. In the proton temperature panel, the red line shows 
the expected temperature following Lopez [1987]. The clock 
angle, 9, is the polar angle in the GSM YZ plane. The 
Alfven speed takes into consideration the influence of the 
alpha particles. The Kan-Lee electric field is defined as: 
Ekl — V X Bt sin^{9j2), where Bt = {By -|- the 

held components being measured in GSM coordinates, and 
9 is the clock angle [see also Shepherd, 2007]. 

There is a CME identihed by Richardson and Cane 
[2010a] starting around 13:00 UT on August 5, 1998 and 
lasting until 12 UT, on August 6. It is a relatively slow 
structure, with an average speed of ~ 380 kms“^ and it 
does not drive a shock ahead of it. The CME does not 
ht the typical characteristics of a magnetic cloud [MC, see 
Burlaga et al., 1981] but it is a relatively clear magnetic 
ejecta. It is given a subjective “mark” of 1 in Richardson 
and Cane [2010a], where 2 represents a clear MC, and 0 a 
very weak and irregular ejecta. It is characterized by a lower- 
than-expected proton temperature, enhanced magnetic held 
strength and relatively smooth rotation of the magnetic held 
vector, among other properties. It has however a relatively 
large density (~ 10 cm“®), and therefore, the proton /?, 
while lower than 1 is only lower than 0.1 in a small part 
of the event. At 07:16 UT on August 6, Wind observes a 
fast forward shock, marked by a red vertical guideline on 
Figure 1, propagating inside this CME. The end boundary 
of the CME is relatively clear around 12:00 UT on August 
6 and corresponds to a short interval of very low magnetic 
held strength and enhanced temperature and P, very similar 
to the interaction region between two CMEs as described in 
Wang et al. [2003a]. Comparing the magnetic held orienta¬ 
tions upstream and downstream of the shock, there is almost 
no change in the clock angle. Overall, these provide strong 
indications that the shock is, indeed, propagating through 
the CME and is not simply at the back boundary, which oc¬ 
curs hours behind the shock. The shock occurs about 18.5 
hours into the CME and about 5 hours before its end. The 
source of the shock is unclear: the rising speed on August 7 
indicates that a fast solar wind speed is catching up to the 
CME. However, the time period around 12 UT on August 
6 is highly reminiscent of the CME-CME interaction region 
as described by Wang et al. [2003a]. 

The conditions upstream of the shock are Np ~ 11 cm“®, 
Vp ~ 370 kms“^, B ~ 10 nT and Tp ~ 10'* K, which cor¬ 
respond to Va ~ 70 kms“* and /? ~ 0.07. The shock is 
quasi-perpendicular, with an angle 0 b„ = 80° ± 5°. The 
density and magnetic field compression ratio are 1.8 ±0.2, 
the shock speed is 460 ±40 kms“*, which corresponds to a 
fast Mach number of 1.5 ±0.3. The Bz component of the 
magnetic field decreases from —6 uT to —11 nT through 
the shock and reaches a peak at —20 nT about 1.3 hours 
behind the shock. The increase in velocity and magnetic 
held strength behind the shock results in a near tripling of 
the coupling function to Fkl = 10 mV m~*, values charac¬ 
teristic of strong driving of the magnetosphere by the solar 
wind. The Dst decreases from —25 nT at 07:00 UT (31 h) 
on August 6 before the shock arrival to —138 nT at 12:00 
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UT (36 h) on the same day. The Sym-H index (equiva¬ 
lent to Dst but with a 1-minute resolution), which is what 
is plotted in Figure 1, drops from —16 uT at 07:15 UT on 
August 6 to —169 nT at 12:14 UT on the same day. This in¬ 
tense geomagnetic storm is quite clearly associated with the 
compression of the magnetic held and the increased veloc¬ 
ity behind the shock. Indeed, the period between the shock 
and the end of the CME is the only time when the cou¬ 
pling electric held, Ekl, reaches high values (>5 mVm“^). 
So the major geomagnetic ehects are related to this shock. 
In fact, the Sym-H decrease, the enhancement of magneto¬ 
sphere convection and two substorms (not shown here) all 
occur in the time period from 8 to 13 UT on August 6. 

This example is representative of our sample of 49 shocks 
inside a CME, as it occurs about 20 hours after the start 
of the magnetic ejecta and propagates into a medium with 
Alfven speed of the order of 80 kms“^ and a plasma /3 of the 
order of 0.1, conditions representative of a relatively weak 
CME. Next, we describe some of the more extreme cases in 
our sample. 

2.2. Very low Mach number shock: November 8, 
1998 

Figure 2 shows Wind measurements around the passage 
at Earth of a shock inside a CME on November 8, 1998. This 
Figure follows the same format as Figure 1. A shock passed 
Wind at 08:15 UT on November 7 followed by a MC-like 
ejecta starting around 22:00 UT (between the black vertical 
guidelines). Starting at that time, the temperature is lower 
than the expected one and the proton /3 is ~ 0.1, which 
are conditions typical of the passage over the spacecraft of 
a magnetic ejecta. At 04:51 UT on November 8, there is 
a fast forward shock propagating through this CME. The 
upstream conditions in the 10 minutes preceding the shock 
passage are: Np ~ 4.5 cm“®, Vp ~ 450 kms“^, B ~ 16.5 nT 
and Tp ~ 4 X lO'* K, which correspond to 14 ~ 170 kms“^ 
and 13 0.03. The shock occurs about 6 hours after the start 

of the magnetic ejecta (as defined by Richardson and Cane 
[2010a]). It is oblique, with an angle Qbv, = 65° ± 10°. The 
density and magnetic field compression ratios are 2 ±0.2, 
the shock speed is 600 ±50 kms“^, which corresponds to a 
fast magnetosonic Mach number of 1.4 ±0.3. The com¬ 
ponent of the magnetic held decreases to —25 nT after the 
shock passage but remains enhanced as compared to the 
upstream for less than 2 hours. The main change in the 
magnetic held is a large increase in the By component. The 
Sym-H index decreases after the arrival of the shock to reach 
two minima 45 minutes and 1.5 hours after the shock pas¬ 
sage with minimum Sym-H values of —173 and —180 nT. 
This corresponds to a decrease of Sym-H by about 40 nT 
from its pre-shock value. Sym-H was decreasing before the 
shock passage and it is impossible to quantify the effects of 
the shock alone. However, the second dip closely following 
the hrst dip is certainly due to the shock effects. 

Note that, even for such a low Mach number, the shock 
compresses the magnetic field and density by a factor of 
~ 2 and may signihcantly affect the geo-effectiveness of the 
CME. As for the previous example, the magnetic held clock 
angle does not vary much through the shock. The enhance¬ 
ment of density downstream of the shock is very short-lived, 
which is a common occurrence for these shocks with very low 
Mach number (it happens for about 25% of the shocks in our 
sample). However, the increases in the magnetic Held and 
pressure (or temperature) are long-lasting. This shock has 
the second strongest upstream magnetic Held and is the one 
occurring closest to the CME leading edge (as a percentage 
of the total CME duration). 

A further example with low Mach numbers is the January 
1, 2006 shock which propagates into the following condi¬ 
tions: Np ~ 5 cm“®, Vp 440 kms“^, B ~ 10 nT and Tp ~ 
510“* K, which correspond to 14 ~ 90 kms“^ and /3 ~ 0.08. 
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Figure 2. Wind observations of the shock inside a CME 
on 11/07/1998. The CME starts around 22:00UT on 
11/07 and is preceded by a shock at 08:15UT on the 
same day (marked with the Hrst dashed red line). A sec¬ 
ond shock passes at 04:51UT on 11/08 (28.85 h), and is 
marked by the second red dashed line. 


The shock is quasi-perpendicular (0s„ = 85° ±5°) with a 
speed of 450 ±50 kms“^, which gives a Mach number of 1.2 
±0.3. The compression ratio for the density and magnetic 
Held are 1.6 ±0.2. As the shock occurs during a northward- 
Bz period, there is no associated geo-effectiveness. 

2.3. Fast Upstream and Alfven Velocities, low 
plasma ( 3 : September 12, 2005 

We next present a case with the largest upstream solar 
wind velocity V = 740 kms“^, which also corresponds to one 
of the largest upstream Alfven speeds. Figure 3 shows the 
plasma and magnetic Held measurements by Wind. There is 
a shock on 01:14 UT on September 11, 2005 which precedes 
a fast CME. The CME duration, as reported by Richardson 
and Cane [2010a] is 26 hours (from 05 UT, September 11 
to 07 UT, September 12), and it is characterized by high 
magnetic Held strength, very low density, low plasma (3, de¬ 
creasing speed proHle. The magnetic Held components are 
not particularly smooth, and this event is marked as 0 in 
Richardson and Cane [2010a]. At 06:00 on September 12, 
Wind encounters a fast-mode forward shock. It occurs at 
or close to the back boundary of the CME, and 25 hours 
after its start. At that time, the magnetic Held is mostly in 
the radial direction, and in addition to the large speed and 
very low density, this leads us to consider that the shock is 
propagating through the rarefaction region behind a CME. 
Development of such rarefaction regions is often seen in nu¬ 
merical simulations [e.g., see Manchester et ai, 2014] and 
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Figure 3. Wind observations of the shock inside a CME 
on 09/12/2005. The CME starts around 05:00UT on 
09/11 and is preceded by a shock at 01:14UT on the 
same day. A second shock passes at 06:OOUT on 09/12 
(30h). 


can also be inferred from remote heliospheric observations 
[Howard et al., 2013]. 

The conditions upstream of the shock are Np ^ 0.5 cm“®, 
Vp ~ 740 kms“^, B ~ 6.5 nT and Tp ^ 9 x 10"* K, which 
correspond to Va 200 kms“^ and P ~ 0.04. The shock 
is quasi-parallel, with an angle 0Bra = 20° ± 5°. The den¬ 
sity and magnetic field compression ratio are 2 ±0.2, the 
shock speed is 1000 ±50 kms“^, which corresponds to a 
fast Mach number of 1.5 ±0.3. The Bz component of the 
magnetic field decreases from —2 nT to —10 nT shortly af¬ 
ter the shock but does not remain southward, consequently 
there are little changes in the Dst index after the shock. 

3. Upstream Conditions and Properties of 
Shocks Inside CMEs 
3.1. Data Selection and List 

We started from the list of CMEs of Richardson and Cane 
[2010a] measured by Wind and compared it with the list 
of shocks measured by Wind as identified and analyzed at 
CfA and the list of shocks measured by ACE as identified 
at UNH by researchers and graduate students there. We 
identified 59 shocks propagating inside or at the back of 
CMEs. In situ measurements were visually analyzed for 
each of these cases, which made us remove 10 cases for which 
the upstream conditions were close to typical, or for which 
the perturbation was not a shock. We ended up with 49 
shocks propagating through CMEs. Our list contains all 


events identified by RC2010 to be shocks propagating in¬ 
side a CME or a CME trailing edge except the following six 
exceptions: one shock propagating inside a potential CME 
(June 14, 2005 marked as CME?), for which we do not con¬ 
sider the upstream medium to be a CME, two shock propa¬ 
gating through CME trailing edges (February 11, 2000 and 
March 31, 2001) for which the upstream conditions are close 
to that of the typical solar wind, and we therefore removed 
from the sample, and three shocks identified by RC2010 as 
propagating inside CMEs (May 4, 1998, May 23, 2000 and 
October 25, 2001). For the first two of these three shocks, 
it was not clear that they satisfy the Rankine-Hugoniot re¬ 
lations for fast-mode forward shocks. The last case had up¬ 
stream conditions comparable to the solar wind, and we also 
removed it from our sample. 

Table 1 lists the 49 shocks propagating through CMEs, 
with the columns showing, from left to right, the CME start 
and end times and type as listed by Richardson and Cane 
[2010a], the shock time at Wind, the delay since the start of 
the CME in hrs, the shock speed in the Earth’s rest frame, 
the shock angle, density and magnetic field compression ra¬ 
tios, the shock Alfvenic Mach number, the upstream fast 
magnetosonic and Alfven speeds, proton /3 and solar wind 
speed, the minimum Dst reached within 12 hours of the 
shock detection at Wind, the drop in Dst from the value be¬ 
fore the shock arrival to the minimum in the next 12 hours 
and the time when the minimum Dst was reached. Out 
of the 49 shocks in our list, 25 were identified in RC2010: 
21 were identified to propagate inside a CME, 1 inside a 
possible CME (marked as CME?), 2 inside a CME trail¬ 
ing edge and 1 inside a CME sheath. We also identihed 24 
additional events; some are clearly not quasi-perpendicular 
shocks, whereas others appear to be and it is unclear why 
they were not in the list of RC2010. There were no shocks 
inside CMEs in 2007 and 2008. For the Bastille day events, 
we used the values provided by Lepping et al. [2001]. Missing 
values are marked with a _. 

Taking into account that two CMEs have multiple shocks 
propagating inside them, and comparing with the full list of 
CMEs measured at 1 AU from Richardson and Cane [2010a], 
we find that about one out of every 7 CMEs has a shock 
propagating through it (47/322). There is no broadly used, 
dedicated shock list; Wang et al. [2010], for example, list 
257 well-defined shocks from 1998 to 2008. Using this value 
for the total number of shocks, we find that about 19% of 
all shocks propagate inside CMEs. This is a lower propor¬ 
tion than that reported by RC2010 for quasi-perpendicular 
shocks (~ 30%). Further discussion is presented in section 
5 about these points. 

In order to compare our dataset with a similar sample 
of shocks propagating into typical solar wind conditions, we 
selected 45 shocks during the years 1997-2006 with approxi¬ 
mately the same yearly distribution as the sample of shocks 
inside CMEs. Hereafter, we refer to this sample of shocks as 
the “control sample”. These shocks were selected from the 
list of well-determined shocks measured by ACE (UNH list) 
and Wind (CfA list). The average properties of this control 
sample are as follows: the fast magnetosonic Mach is 2.34 
(average Alfven Mach of 2.84), the average shock speed is 
540 kms“^, the average p is 0.63 (median 0.49), the aver¬ 
age compression ratio is 2.55 (median 2.4). This sample has 
average and median properties close to but slightly higher 
than that found by other researchers for a larger number of 
shocks. For example. Oh et al. [2007], found a typical com¬ 
pression ratio of 2, based on 249 shocks from 1995 to 2001, an 
Alfven Mach of 2.1, a typical shock speed of 490 kms“^, and 
a upstream P ~ 0.3. From the list of Wang et al. [2010], the 
average speed of the 257 shocks is 500 kms“^, the average 
fast magnetosonic Mach is 2.06 and the average compression 
ratio is 2.06. The larger values for our control sample may 
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Figure 4. Superposed epoch analysis of the 49 shocks with data from omniweb. The left panels show 
the velocity, proton density, proton temperature, total magnetic field, absolute value of the magnetic 
field Bz component in GSM coordinates, Alfven speed and Sym-H index, from top to bottom. The right 
panels show the dimensionless velocity, density, temperature, total and Bz magnetic field, from top to 
bottom. 


be due to the fact that it includes more shocks during the 
rise and solar maximum phases (1999-2002) than a random 
sample would have. Shocks at solar maximum tend to be 
faster and have larger Mach numbers and compression ratio 
[Oh et ai, 2007]. This may also reflect the fact that this sam¬ 
ple excludes shocks inside CMEs which tend to have lower 
compression ratio and Mach numbers, as discussed below. 

3.2. Superposed Epoch Analysis 

We performed a superposed epoch analysis of the 49 
shocks propagating inside a CME, using 1-minute data from 
omniweb^. Omniweb data is time-shifted to the nose of 
Earth bow shock, providing a common temporal reference 
point for all the shocks. For the November 6, 2001 event, 
there was no data after the shock arrival (except Sym-H) 
and for the July 15, 2000, there was no data except Sym-H 
in the considered time window. These were associated with 
the two strongest drops of Sym-H after the shock from our 
dataset, so we included them nonetheless in our analysis. 
The October 31, 2001 has about 3.5 hours of data missing 
around the shock arrival time, and it is also included in the 
sample. We visually identify the shock arrival time to serve 
as the zero epoch. Figure 4 shows the velocity, density, tem¬ 
perature, magnetic field strength, absolute value of GSM Bz 
component, Alfven speed and Sym-H index for 4.5 hours be¬ 
fore and after the shock arrival time. Because a few events 
with extreme values may affect the average of the events, we 
also plot in the right panel of Figure 4, the velocity, density, 
magnetic field strength, GSM Bz component and temper¬ 
ature in dimensionless form. To do this, we calculate, for 
each event, the maximum value of the parameters before 
the shock arrival and we scale the measurements by that 


value before doing the averaging. The averaged results are 
then scaled so that the pre-shock value is about 1. In this 
way, any event with extreme values of the plasma and mag¬ 
netic parameters shall have less impact on the global trend. 
It is also easy to visualize the approximate average jump at 
the shock with this plot. The same plots are shown for the 
control sample in the Appendix. 

The left panel of Figure 4 clearly shows the average up¬ 
stream conditions into which the shocks propagate: velocity 
of 450 kms“^, density of 6 cm“®, magnetic field of 8.5 nT 
(with a Bz component of about 4.5 nT), and an Alfven speed 
of about 105 kms“^. These average values are very close 
to the median (average) from our sample as listed in Ta¬ 
ble 1, which are 430 (455) kms“^ 4.5 (5) cm“^ 8 (8.5) 
nT, a Alfven speed of 85 (105) kms“^ and a plasma /3 of 
0.15 (0.08). These upstream conditions differ from typical 
ones as found for our control sample: upstream velocity of 
410 kms“^, magnetic field of 6 nT and Alfvenic speed of 
70 kms“^. A difference of 30-40 kms“^ in both upstream 
solar wind and Alfvenic speeds may not sound like much, 
but, for a shock with a speed of 600 kms“^, the Alfvenic 
Mach number would be 2.7 for normal conditions but only 
1.4 for upstream conditions as encountered inside CMEs. 

As shown in the right panel of Figure 4, the shock com¬ 
presses, on average, the density by a factor of 2, and the 
magnetic field by a factor of 1.7. The average increase in 
the magnitude of the Bz component is about 2. As shown 
in the left panel of Figure 4, the magnetic field continues to 
increase behind the shock, so that the maximum magnetic 
field strength is about 17.5 nT, 1 hour behind the shock, 
and the average density increases to about 13.5 cm“®. The 
jump through the shock in the velocity is about 100 kms“^. 
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Figure 5. Left panel: yearly distribution of the shocks inside CMEs (blue bars) and shocks from the 
control sample (yellow bars) as well as total number of CMEs measured in situ (red triangles). Right 
panel: Shock speed and Aifven Mach number for the shocks inside CMEs (blue diamonds) and shocks 
from the control sample (yellow circles). 


Although the events were not selected based on their 
geo-effectiveness, the average Sym-H index is perturbed 
at a value of —30 nT before the shock arrival. This is 
somewhat consistent with the fact that the upstream 
conditions correspond to the passage at 1 AU of a CME 
[see for example Richardson and Cane, 2010a]. On av¬ 
erage, there is a sudden storm commencement at the 
shock passage with an increase of the Sym-H by about 
14 nT, followed by a decrease with the Sym-H reaching 
—51 nT 4.5 hours after the shock passage. Note that the 
inclusion of the two strongest drops for which there is no 
downstream plasma data does not change the shape of 
the Sym-H curve and only changes the minimum value 
by 3 nT. The shape of the Sym-H curve is very similar 
to that found for the control sample, but it is shifted by 
—20 nT. Therefore, a moderate geomagnetic storm is 
the average result from the passage at Earth of a shock 
inside a CME. 

3.3. Detailed Analysis 

We now give some additional details regarding the 
upstream conditions and the geo-effectiveness of these 
shocks. The two subsequent sections discuss the shock 
properties and location within the CME in more detail. 

20 of the 49 shocks propagate inside solar wind 
with Alfvenic speed greater than 100 kms“^. Some of 
the largest Alfvenic speeds are measured for the three 
shocks during the Halloween 2003 time period (October 
26-28) with upstream Alfvenic speeds between 180 and 


220 kms“^. The March 20, 2002 shock has upstream 
Alfvenic speed in excess of 400 kms“^ due to very low 
upstream density. In fact, 19 shocks propagate through 
upstream densities lower than 3 cm“^, 13 of which also 
have Alfvenic speed greater than 100 kms“^. 

Only 15 shocks propagate inside magnetic fields 
stronger than 10 nT, which is the average magnetic 
field measured at 1 AU inside CMEs [Richardson and 
Cane, 2010a]. 3 of these propagate inside magnetic 

fields stronger than 15 nT, with the November 6, 2001 
shock the one propagating into the strongest magnetic 
field (~ 20 nT, although the upstream density is greater 
than 10 cm In that sense, many of the CMEs into 
which the shocks propagate are relatively weak. Most 
of the CMEs are also relatively slow. 8 shocks occur 
with upstream speeds greater than 550 kms“^ with the 
largest upstream speed of 740 kms“^ as already dis¬ 
cussed for the September 12, 2005 shock. It is expected 
that shocks are neither able to overtake extremely fast 
CMEs, nor able to propagate through CMEs with such 
fast upstream conditions. This might explain why the 
median upstream speed (430 kms“^) for the 49 shocks 
is slightly lower than the median speed of CMEs mea¬ 
sured at I AU [476 kms“^, see: Richardson and Cane, 
2010a]. 

The average (median) Dst measured within 12 hours 
of the shock arrival at Wind is —81 (—72) nT. 19 of 
the 49 shocks are associated with an intense geomag¬ 
netic storm (Dstpeak < —100 nT) within 12 hours of 
the shock detection at Wind. The two largest storms 
are the one on July 15, 2000 where Dst drops from 
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Figure 6. Compression ratio of the shocks vs Alfven Mach number (left panel) and plasma (3 (right 
panel) for the shocks inside CMEs (blue diamonds) and from the control sample (yellow circles). 


—55 nT to —289 nT in the 8 hours following the shock 
arrival and the one on November 6, 2001 where the 
Dst drops from —73 nT to —292 nT 5.5 hours after the 
shock arrival. 15 shocks were associated with a drop of 
the Dst index of more than —50 nT within 9 hours of 
the shock passage. The average drop within 12 hours is 
—46 nT (median —32 nT). Since these shocks occur in¬ 
side CMEs, it is not always possible to separate exactly 
the effects of the CME with that of the compression 
from the shock. However, considering the typical jump 
at the shock, the dawn-to-dusk electric field increases on 
average by a factor of 2.5 through the shock. The large 
drop in Dst found within 9 hours of the shock passage 
for 30% of the shocks is also consistent with the com¬ 
pression due to the shock (and possibly the following 
CME or fast solar wind stream) to contribute signifi¬ 
cantly to the additional geo-effect. In short, for some 
events, such as the one shown in Figure 1, there are ge¬ 
omagnetic effects clearly caused by the shock passage, 
and they can be very intense. 

3.4. Shock Properties 

The left panel of Figure 5 shows the yearly distribu¬ 
tion of the shocks as compared to that of CMEs. The 
number of shocks itself is known to have a solar cycle 
dependence [Oh et al., 2007]. This effect is probably 
amplified for shocks inside CMEs since the number of 
CMEs also has a solar cycle dependence. Nonetheless, 
we found 2 occurrences before mid-1998 and 8 after mid- 
2004, which are time periods of moderate activity when 
the monthly sunspot number was less than 50. 

The right panel Figure 5 shows the distribution of 
Alfven Mach numbers and shock speeds for the shock 
inside CMEs and for those from he control sample. Typ¬ 
ical error for the Mach number and the compression ra¬ 
tio is ±0.3, it is ±70 kms“^ for the shock speed. From 
this Figure, it is clear that many of the shocks prop¬ 
agating through CMEs have a very low Alfven Mach 
number (45% have Ma < 1.75 compared to 7% for 
our control sample). However, very few of them have 


slow speeds (only 1 shock propagating through a CME 
has a speed of less than 450 kms“^, as compared to 
9 from the control sample). The average (resp. me¬ 
dian) speed of the shocks propagating through CME is 
617 (resp. 580) kms“^, significantly larger than that 
from the control sample (average speed of 540 kms“^), 
and much larger than that from the other studies cited 
above. The average (resp. median) Alfven Mach num¬ 
ber is 2.3 (resp. 1.9), which is slightly lower than that 
from other samples (average of 2.9 for the control sam¬ 
ple). In summary, these 49 shocks propagating inside 
CMEs are on average fast but weak, which is what was 
expected. This comes from the following facts: (i) only 
fast shocks can overtake a previous CME and propagate 
through it, and (ii) as the upstream and Alfven speeds 
are high inside the CME, these shocks are weak. As 
noted in RC2010, what is referred to as the shock prop¬ 
erties (shock angle, Mach numbers) is, in fact, strongly 
dependent on the upstream conditions. 

Figure 6 shows the compression ratio with respect 
to the Alfven Mach number and the compression ratio 
with respect to /3. The diamonds show the 49 shocks 
propagating through CMEs, whereas the circles show 
the 45 shocks from the control sample. It appears that, 
at low Aflven Mach numbers, shocks propagating inside 
CMEs have a slightly larger compression ratio as com¬ 
pared to shocks from the control sample at the same 
Mach number. We come back to these considerations 
in the section 4 to show how the compression ratio of a 
shock increases with decreasing upstream /3. 

3.5. Shock Location 

In the list of CMEs of Richardson and Cane [2010b], 
the median duration of a CME is about 30 hours. As 
shown in the left panel of Figure 7, only about one 
third of the shocks occur within 16 hours of the start 
of the CME. In fact, one third occurs more than 32 
hours after the CME start. We also found no case for 
which the shock propagates within 3 hours of the start 
of the CME. The extreme cases correspond to very large 
CMEs, as for example the April 17-19, 2002 CME (du- 
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Figure 7. Distribution of the shock position inside the CME. In the left panel, the delay between the 
start of the CME and the shock location is plotted in hours; the right panel shows the location of the 
shock as a proportion of the CME total duration, where 1 indicates a shock at or after the back boundary. 


ration 47 hours, shock after 40.5 hours) and the Oc¬ 
tober 3-5, 2000 (duration 41 hours, shock at the back 
boundary). The three most extreme cases probably cor¬ 
respond to shocks propagating through complex ejecta 
composed of multiple CMEs: the August 20-23, 1998 
CME (duration of 60 hours with a shock after 48 hours 
and a shock at the back boundary) and October 22-24, 
2003 CME (duration of 61 hours with a shock at the 
back boundary); both have a radial size of more than 
0.7 AU or twice more than a typical CME. 

The right panel of Figure 7 shows the location of 
the shock within the CME as a proportion of the CME 
duration, where 0 identifies a shock at the CME front 
boundary and 1 a shock at (or after) the CME back 
boundary. 38% (19/49) of the shocks occur at the back 
boundary of the CME. Even removing these, the dis¬ 
tribution is clearly biased towards shocks being close 
to the back of the CME: the median (average) position 
of the 30 other shocks is 80% (72%) from the CME 
front with the closest shock 25% of the way from the 
CME front (Nov. 8, 1998, see Figure 2). Only 7 shocks 
are in the front half of the CME, and 11 in the front 
70%. It should be kept in mind that the part of the 
CME downstream of the shock has been compressed, 
as shown for example in Lugaz et al. [2013]. Therefore, 
a shock occurring temporally halfway through the CME 
has already propagated through more than half of that 
CME. However, considering that the compression ratio 
of the shock is typically low, and that so few shocks are 
found even in the front 70% of the CME, we believe 
that the distribution is truly biased. 

It is our opinion that this biased distribution reflects 
the fact that shocks weaken due to the large Alfvenic 
and upstream speeds inside CMEs and only some shocks 
“survive” the entire propagation throughout a CME. 
Due to CME expansion, the upstream medium is typ¬ 
ically faster by 100-200 kms“^ at the front boundary 
of the CME as compared to the back. Many of the 
shocks in our sample would not be fast-mode shocks 
if the upstream speed increased by 100 kms“^, every¬ 
thing else being equal. Looking at Figure 3, it is very 
unlikely that the September 12, 2005 shock would re¬ 


main a fast-mode shock in the front half of the CME 
where the upstream speed is greater than 1,000 kms“^, 
i.e. greater than the shock speed. Another example is 
shown in Figure 8 with the September 30, 2001 shock. 
There, a shock propagates at the back of a fast and large 
CME, about 43 hours after the CME start and 5 hours 
before its end. The shock speed at 700 kms“^ is lower 
than the speed of the CME front (^ 750 kms“^). Addi- 



Time (hr) sfnc® 09/2^ OOUT 


Figure 8. Wind observations of the shock inside a CME 
on 09/30/2001. The CME starts around lliOOUT on 
09/29 and is preceded by a shock at 09:40UT on the 
same day. A second shock passes at 19:14UT on 09/30 
(43.2 h). 
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tionally, the Alfven speed also increases by 100 kins“^ 
to 150 kms“^ from the back to the center of the CME. 
With a shock speed of 700 kms“^, it is clear that the 
shock will not remain a fast-mode shock as it reaches 
the center of the CME. An example of a fast-mode wave, 
which probably is the result of a fast-mode shock decay¬ 
ing as it propagates through a CME was given in Wang 
et al. [2003b] for the October 03-05, 2001 CME (their 
Figure 9). Counter-examples are the shocks shown in 
Figures 1 and 2 , which are likely to remain fast-mode 
shocks throughout the CME. 

We should also note that there are some measure¬ 
ments of shocks propagating through a CME sheath (for 
example the March 30, 2001 sequence of two shocks), 
which may indicate that some shocks can indeed prop¬ 
agate through an entire magnetic ejecta. Overall, our 
sample suggests that not all shocks can remain fast¬ 
mode shocks as they propagate through a CME. A 
rough estimate of the number of “missing shocks” can 
be obtained as follows. First, we cannot assume that 
shocks are as likely to be measured in all parts of CMEs. 
This is simply because, when a shock is measured in¬ 
side a CME, the back portion of this CME (downstream 
region of the shock) has gone through compression, 
which shortens its duration. Consequently, it changes 
the probability of detecting a shock at a given location 
as measured from the CME total duration. For exam¬ 
ple, the shock on May 29, 2005 is measured half-way 
through a CME 6 hours after its start and 6 hours be¬ 
fore its end. However, the 6 hours behind the shock, 
correspond to compressed CME material, not the “ac¬ 
tual CME” before compression. This means that the 
shock has already propagated through more than half 
of the CME towards the front. 

The typical compression ratio, in density and mag¬ 
netic field, of the shocks in our study, is 2. Therefore, 
we assume that the duration of the compressed portion 
of these CMEs has been reduced by a factor of 2 as com¬ 
pared to the duration before the shock passage. Then, 
under this assumption, we can estimate that an equal 
number of shocks should be measured in the front 67% 
of the CMEs as in the back 33% (based on the CME 
duration). 11 shocks were found in the front 67%, 19 
shocks at the CME back boundary and 19 shocks in 
the back 33% (without including the shocks at the back 
boundary). We estimate that 25±10 shocks should have 
been measured in the front 67% of the CME for the 
same number of shocks measured at the back of the 
CME. We can therefore give a rough estimate of 50% 
± 25% of shocks which do not remain fast-mode shocks 
during propagation through an entire CME. 

4. Analytical Considerations 

In this section, we discuss the compression of a MHD 
shock under different upstream conditions using the an¬ 
alytical Rankine-Hugoniot formalism. Most of this work 
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Figure 9. Compression ratio (in density) of a fast-mode 
shock assuming the shock normal is anti-parallel to the 
velocity. The first two panels show the cases of a perpen¬ 
dicular shock, while the last panel shows the dependence 
of the compression ratio on the shock angle, Alfven Mach 
number and plasma p. 

will be well known to researchers working on Earth’s 
bow shock but these equations have been rarely used 
to discuss CME-driven shocks. For the sake of simplic¬ 
ity, we limit ourselves in the next section to shocks for 
which the normal is anti-parallel to the velocity vec¬ 
tor. As CME-driven shocks tend to be relatively blunt 
[see for example Janvier et al., 2014], most shocks are 
observed with a normal within 30° of the GSE x di¬ 
rection, which is, in general, the direction of the speed. 
This assumption is therefore justified for these simple 
analytical considerations. 

It is important to remember that the plasma P is re¬ 
lated to the ratio of the Alfven, Ma to sonic. Mg Mach 
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numbers as follows: 

M2 2 


( 1 ) 


Under this assumption, the formula for the compression 
ratio, \/X of a MHD shock is the positive root of the 
cubic equation [Zhuang and Russell, 1981; Cairns and 
Grabbe, 1994]: 


aX^ - bX^ + cX - d = 0, (2) 

with: 

a = (7 + 1)M1 

b = (7-1)M^ + (7 + 2)M1cos2 6» + 7(1 + /3)M1 
c = (7 — 2 + 7 cos^ 6 *)M ]4 + (7 + 27/3 + 1)M^ cos^ 0 
d = (7 — 1 )M^ cos^ 0 + 7 / 3 cos"^ 0 

where Ma is the Alfven Mach number and 0 = 0Bn the 
angle between the shock normal and the upstream mag¬ 
netic field vector. The full version of this equation is a 
more complex cubic equation as discussed in Petrinec 
and Russell [1997] and Kabin [2001], which involves the 
angle between the shock normal and the velocity vector. 

To start the discussion, we can make a further simpli¬ 
fication, assuming a perpendicular shock (^Bn = 90°). 
The equation simplifies to a quadratic equation, since 
d = 0. For 7 < 2, the equation always has a unique 
positive root, X±. However, the shock must also be 
compressive (Xj_ < 1 ) to be an acceptable solution to 
the Rankine-Hugoniot relations, which can be written 
as: 


(7 + l)Ml - (7 - l)Ml - q(l + /3) + 7 - 2 > 0 


or 

Ml > 1 + ^ (3) 

Equation 3 clearly shows that MHD shocks with ex¬ 
tremely low Alfven Mach numbers can only be observed 
when the plasma /3 is low. For example, for 7 = 1.5, 
the minimum Alfven Mach number at /3 = 1 is 1.32; 
whereas for /3 = 0.1, it is 1.04. Below this limit, the 
shock fast magnetosonic Mach number becomes less 
than 1. Talking only about the fast magnetosonic Mach 
number or the Alfven Mach number hides the fact that 
some of these shocks have large sonic Mach numbers. 
The left and middle panels of Figure 9 show the com¬ 
pression ratio of a perpendicular shock for different val¬ 
ues of /3 and the Alfven Mach number, illustrating how 
weak shocks may be significantly more compressive for 
/3 < 0.1. 

Examples of the solutions to equation (2) for oblique 
shocks are shown in the last panel of Figure 9, where 
the compression ratio is shown for different values of 
/3 = [0.1, 0.4,1.2] and for Ma = [1.5, 2, 3]. For /3 = 0.01, 
the values are nearly identical to that for jS = 0.1. Note 
also that the sonic Mach number is the same for the 


solid black curve (/3 = 0.1, M^ = 1.5) and the dotted 
red curve (/3 = 0.4, M^ = 3), and the fast Mach num¬ 
ber is about the same between the solid black curve 
and the blue dashed curve (/3 = 1.2, Ma = 2). In a 
low-/3 plasma, a shock with an Alfven Mach of 1.5 has 
a compression ratio between 1.6 and 2 , depending on 
the shock angle. For /3 = 1, the compression ratio is 
between 1.1 and 1.6. This is somewhat consistent with 
the behaviors found for the 49 shocks inside CMEs as 
compared to the 45 shocks in more typical upstream 
conditions, as shown in the left panel of Figure 6 . 

The compression of the tangential component of 
the magnetic field is given by l/Aut with Aut = 
X clear that the tangential magnetic 

compression is always larger than the compression in 
density, the two being equal for perpendicular shocks. 
For small Alfven Mach and quasi-parallel situation, the 
solution tends towards the switch-on shock. It should 
be remembered that for shocks inside CMEs, the shock 
angle depends primarily on the local orientation of the 
magnetic field. It is likely that the shock angle dramat¬ 
ically changes as the shock propagates inside the CME. 

5. Discussions and Conclusion 

In this article, we identified 49 shocks propagating 
inside or closely following a CME during solar cycle 23 
(1997-2006, with no such shock during 2007-2008). We 
analyzed the properties of these shocks and compared 
them to 45 shocks propagating into normal solar wind 
conditions. We found that shocks propagating inside 
CMEs are found to be, as expected, fast and weak with 
a median speed of 580 kms“^ and a median Alfven 
Mach number of 1.9. These shocks represent about 15- 
20% of all shocks measured at 1 AU during solar cycle 
23. RC2010 reported that 30% of quasi-perpendicular 
shocks during the same time period were found inside or 
following CMEs. About 70% of the shocks in our sample 
have a shock angle between 65° and 90°, i.e. are quasi¬ 
perpendicular (however, 24 out the 49 shocks analyzed 
here were not part of the RC2010 study). As detailed 
in RC2010, the fact that shocks inside CMEs repre¬ 
sent a larger proportion of quasi-perpendicular shocks 
than that of the whole shock population is a direct con¬ 
sequence of the conditions required for a shock to be 
quasi-perpendicular at 1 AU. If a spacecraft crosses a 
shock close to its nose, quasi-perpendicular conditions 
can only be met if the magnetic field is primarily in the 
y — z plane. This is often the case inside a CME but 
rarely outside of one. When the “wing” of a shock is 
measured, the shock may be quasi-perpendicular under 
normal solar wind conditions. However, such measure¬ 
ments are rarer [Janvier et at, 2014]. In short, shocks 
inside CMEs have a higher likelihood of being quasi¬ 
perpendicular than random shocks. 

Compared to the full list of CMEs measured at 1 AU, 
we find that about 9% of CMEs have a shock propagat¬ 
ing through them, and an additional 6 % have a shock at 
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the location of the back boundary [following the bound¬ 
aries selected by Richardson and Cane, 2010a]. This 
proportion varies significantly with solar cycle. For 
2003, when the total number of CMEs measured at 
Earth was relatively low, but many occurred during 
a short time period (the October-November months, 
the so-called Halloween storms), as many as 30% of the 
CMEs identified at Earth had a shock inside them. 

Large southward magnetic fields are the leading cause 
of geomagnetic storms, and these occur preferentially 
inside CMEs. Here, we found that 19 occurrences of 
shocks inside CMEs were associated with an intense 
geomagnetic storm within 12 hours of the shock pas¬ 
sage. This represents over 20% of all intense geomag¬ 
netic storms identified in solar cycle 23. Shocks inside 
CMEs are therefore an important cause of intense geo¬ 
magnetic storms at 1 AU. In four cases, the downstream 
magnetic field a few hours after the shock was larger 
than 30 nT. In addition, the August 24, 2005 CME has 
a sharp discontinuity 3 hours after the shock, which 
compressed the magnetic field from 25 nT to close to 
60 nT. These large values of the magnetic field strengths 
are comparable to that of CMEs driving superstorms 
(Dstpeak < —250 nT). However, here they are caused by 
weak-to-average CMEs. This shows that shocks propa¬ 
gating inside CMEs are one of the ways to create strong 
Bz (and strong magnetic Held) by compressing weak or 
average CMEs into CMEs with more extreme values of 
the magnetic field strength. 

We did not systematically investigate the source of 
the shocks propagating inside CMEs, but we believe 
that most of them are CME-driven shocks for the stud¬ 
ied period. First, most shocks at 1 AU are CME-driven. 
Eor example, Jian et al. [2006a] and Jian et al. [2006b] 
reported 151 CME-driven shocks from 1995 to 2004 
as compared to a total of 67 forward shocks driven 
by SIRs or CIRs during the same time period. Oh 
et al. [2007] found a lower proportion of CIR-driven 
shocks at about 16% from 1995 to 2001. They also 
found that CIR-driven shocks at 1 AU were signifi¬ 
cantly weaker and slower than CME-driven shocks (420 
as compared to 500 kms“^). This makes it unlikely 
that CIR-driven shocks are able to overtake and prop¬ 
agate through CMEs, although they may contribute to 
the large number of shocks propagating at the back of 
CMEs. In addition, we find that many CMEs into which 
the shocks propagate are complex (for example, about a 
quarter of them have a duration of more than 30 hours), 
which may indicate that they are in fact the result of 
the interaction of multiple CMEs. 

We have also analyzed the location of the shocks 
within the CMEs. Most shocks occur more than 20 
hours after the start of the CME in the back half of 
the magnetic ejecta. Even taking into consideration the 
compression of the back of the CME by the shock, it 
appears that fewer shocks are measured in the front 
part of the CME as what could be expected from an 


uniform distribution. We conclude that about half of 
the shocks entering a CME may not remain fast-mode 
shocks throughout the CME. This is likely due to the 
larger upstream solar wind speed, due to the CME 
expansion, and the location of the maximum in the 
Alfvenic speed, which often occurs close to the mid¬ 
dle of the CME. In the CME frame, the shock is often 
near-stationary, while the back of the CME propagates 
sunwards towards the shock due to the CME expan¬ 
sion. Future work should further investigate possible 
signatures of shocks after they exit CMEs. 
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Appendix A: Control Sample 

The yearly distribution of the 45 shocks from our con¬ 
trol sample is shown in the left panel of Figure 5. The 


main properties of these shocks are shown in Figures 5 
and 6. Figure 10 shows the superposed epoch analy¬ 
sis for these 45 shocks in the same format as Figure 4. 
Table 2 lists the properties of these shocks. 
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Figure 10. Superposed epoch analysis of the 45 shocks from our control sample. The left panels show 
the velocity, proton density, proton temperatnre, total magnetic field, absolnte valne of the magnetic 
field Bz component in GSM coordinates, Alfven speed and Sym-H index, from top to bottom. The right 
panels show the dimensionless velocity, density, temperatnre, total and Bz magnetic field, from top to 
bottom. 
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Notes 

1. http://oiniiiweb.gsfc.nasa.gov 
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2.2 

1.95 

112 

107 
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Median 



22 
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2.1 

1.9 
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Table 1: List of 49 shocks propagating inside CMEs from 1997 to 
2006. The columns list the CME start and end times, and type as 
given by Richardson and Cane [2010a] (2 = MC, 1 = MC-like, 0 
= irregular CME); the shock time, the time delay (in hours) since 
the start of the CME (bolded At indicates that the shock time 
occurred before the CME end time), the shock speed, Alfven Mach 
number, angle, compression ratios in density and magnetic field; 
the upstream fast magnetosonic, Alfven speed, f3 and solar wind 
speed; the minimum Dst and drop in Dst in the 12 hours following 
the shock passage and the delay in hour between the shock arrival 
and the minimum Dst. 
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480 

2.3 

15 

4 

2.4 

65 

0.15 

340 

08/29 10:16 

480 

2.1 

30 

1.9 

1.5 

40 

0.5 

400 

11^7 18 :28 

750 

3 

60 

2.3 

2.3 

100 

450 

1 

2005 

06/14 18:35 

540 

4.5 

65 

2.1 

1.8 

30 

0.8 

435 

07/10 03:44 

520 

2 

80 

1.8 

1.8 

90 

0.1 

350 

08/01 06:47 

490 

1.7 

85 

1.9 

1.9 

95 

0.15 

500 

09/02 14:24 

580 

2.4 

65 

2.1 

2.1 

60 

0.6 

425 

09/15 08:52 

660 

2.4 

55 

2.8 

2.5 

60 

0.4 

535 

2006 




LUGAZ ET AL.: SHOCKS INSIDE CMES: 1997-2006 


X - 17 


09/04 00:12 505 2.9 45 


2.2 1.7 


45 


0.7 405 


Table 2: List of 45 shocks from our control sample from 1997 to 
2006. The columns list the shock time, shock speed, Alfven Mach 
number, angle, compression ratios in density and magnetic held; 
the upstream Alfven speed, /3 and solar wind speed. 



